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were  6 in.  in  diameter  with  a 6-in.  long,  ogive  nose  (CRII  = 1.25).  Two 
penetrator  weights  were  selected,  141.4  Ib  (Case  1)  and  424.1  lb  (Case  2),  to 
give  projectile  sectional  pressures  of  5 psi  for  Case  1 and  IS  psi  for  Case  2.f^. 

At  impact,  the  force  loadings  on  the  rigid  penetrators  for  the  two 
cases  arc  the  same,  = ^2  ’ ratio  of  decelerations  will  be  inversely 

proportional  to  the  penetrator  weights,  i.e.,  aj/a2  ~ ~ ^ 

penetration  proceeds,  the  velocities  will  differ,  due  to  the  different  de- 
celerations, and  the  ratio  of  forces,  ^2^^!  ’ scaled  deceleration  ratio, 

^2'''2^^1*'*1  ’ depart  from  unity. 

Each  penetrator  problem  was  run  from  initial  impact  to  a penetration 
depth  of  8 in.,  or  2 in.  below  full  nose  embedment.  From  the  solution  results, 
the  scaled  deceleration  ratio,  a.,W.,/ajWj  , was  computed  as  a function  of  depth 

of  penetration.  These  results  showed  that  the  scaled  deceleration  ratio 
oscillated  around  a value  of  1 for  about  the  first  half  of  nose  embedment; 
the  ratio  then  steadily  increased  to  a value  of  1.06  at  a penetration  depth 
of  8 in.  The  gradually  larger  force  applied  to  the  heavier  projectile 
(Case  2)  is  due  to  its  penetration  velocity  staying  higher  (due  to  its  lesser 
deceleration),  relative  to  Case  1. 

Oetailcd  results  of  the  numerical  solutions,  including  time  histories 
of  the  penetration  dynamics  variables,  force  loading  distributions,  and  field 
plots  of  the  target  response,  are  presented  in  the  report. 
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This  report  was  prepared  by  Messrs.  Mark  II.  Wagner  and 
Christopher  C.  Fulton  of  California  Research  and  Technology,  Inc., 
Woodland  Hills,  California, under  PO  DACA  39-77-M-0136  as  a part  of 
ongoing  work  at  the  U.  S.  Army  Engineer  Waterways  Experiment  Station 
(WES),  PA  Project  No.  4A161102AT22,  Task  A2,  Work  Unit  006,  "Effective- 
ness of  Earth  Penetrators  in  Various  Geologic  Environments."  The 
principal  investigator  at  California  Research  and  Technology  was 
Mr.  Wagner.  The  programming  and  computer  runs  were  conducted  by 
Mr.  Fulton. 

The  problems  were  specified  and  the  work  was  monitored  by 
Dr.  Behzad  Rohani  of  the  Soil  Dynamics  Division  (SDD) , Soils  and 
Pavements  Laboratory  (S5PL) , WES,  under  the  supervision  of 
Dr.  J.  G.  Jackson,  Jr.,  and  the  general  direction  of  Mr.  James  P. 

Sale,  Chief  of  SfiPL. 

COL  J.  L.  Cannon,  CE,  was  Commander  and  Director  of  WES 
during  the  period  this  work  was  performed.  Mr.  F.  R.  Brown  was 
Technical  Director. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units 
units  as  follows: 

Multiply 

feet  per  second 
inches 

kips  (force) 

microns 

pounds  (mass) 

pounds  (mass)  per 
cubic  foot 

pounds  (mass)  per 
square  foot 


of  measurement  can  be  converted  to  metric  (SI) 


By 

To  Obtain 

0.3048 

metres  per  second 

25.4 

millimetres 

4.448222 

kilonewtons 

0.001 

millimetres 

0.4535924 

kilograms 

16.01846 

kilograms  per 

cubic  metre 

4.882428 

kilograms  per 

square  metre 
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NUMERICAL  ANALYSES  OF  PENETRATION  DYNANtICS  IN  SUPPORT  OF 
I NXTSTI CATIONS  OF  SCALING  REUTIONS  FOR  E.ARTU  PEMITRATORS 

1 . SCOPE 

This  report  describes  the  results  of  two  f i ni te -di fference 
code  calculations  of  projectile  penetration  conducted  by  CRT 
to  provide  information  of  use  in  the  study  of  penetrator  scaling 
relations.  These  two  calculations  were  designed  specifically 
to  assess  the  effect  of  projectile  sectional  pressure  (W/A) 
on  penetration  dynamics. 

2.  PENETRATOR  PROBLEM- -DESCRIPTION 

The  penetrator  problem  considered  was  the  normal  impact 
of  a rigid  body  projectile  into  a sand  target  at  200  ft/sec, 
as  shown  in  Figure  1.  The  conditions  of  the  two  cases  con- 
sidered are  listed  in  Table  1.  The  projectile  for  Case  2 
weighed  3 times  as  much  as  the  Case  1 projectile,  giving  a 
projectile  sectional  pressure,  W/A,  of  5 psi  for  Case  1 and 
15  psi  for  Case  2;  all  other  conditions  of  the  problems  were 
identical . 

The  nominal  properties  of  the  sand  target  are  listed  in 
Table  2.  The  material  model  used  was  the  same  as  that  used 
for  Layer  2 of  the  Watching  Hill  Site  at  DRES  in  a previous 
study\  except  that  the  low-pressure  end  of  the  yield  surface 
was  modified  so  that  the  cohesion  would  be  equal  to  zero. 

The  penetrator/target  interface  was  assumed  to  be  friction- 
less for  these  problems. 

3.  COMPUTATIONAL  METHOD 

Numerical  solutions  of  the  penetrator  problems  were 
obtained  with  WAVE-L,  a two-dimensional  finite-difference 
Lagrangian  hydrodynamic-elastic-plastic  code.  This  code  has 
been  employed  to  solve  penetrator  problems  on  a number  of  DNA 
and  WES  projects  during  the  past  few  years. 
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The  computational  grid  set  up  for  these  problems  is 
partially  shown  in  Figure  2.  The  basic  cell  size  used  in  the 
target  was  Ar^  = AZq  = 0.5  in.,  or  1/b  the  projectile  radius. 
Beyond  a radius  of  6 in.  and  a depth  of  12  in.,  the  cell 
dimensions  were  gradually  increased,  in  2%  steps.  The  over- 
all grid  extended  to  a radius  of  30.5  in.  and  a depth  of 
36.5  in.,  comprising  2800  grid  points.  This  grid  extends 
beyond  the  boundaries  of  wave  propagation  that  took  place 
during  these  problems  (which  were  run  to  a penetration  depth 
of  8 in.),  so  that  they  could  be  restarted  and  continued  to 
longer  times  if  desired. 

The  integration  time  step  varies  during  the  course  of 
the  problem,  depending  on  the  sound  speed,  minimum  cell 
dimension,  and  stability  criteria  selected.  For  these  problems, 
the  average  time  step  was  1.5  psec. 


4.  NUMERICAL  SOLUTION  RESULTS 


Each  penetrator  problem  wa 
a penetration  depth  of  8 in.,  o 
ment.  The  computed  penetrator 
tration  for  each  case  is  shown 
the  penetrator  deceleration,  ve 
and  attached  length  (vertical  d 
point  of  separation  of  projecti 
Figure  4.  The  deceleration  for 
mately  3 times  that  of  Case  2, 
inversely  proportional  to  the  p 
This  scaling  ratio  can  be  exami 


s run  from  initial  impact  to 
r 2 in.  below  full  nose  embed- 
deceleration  vs  depth  of  pene- 
in  Figure  3.  Time  histories  of 
locity,  depth  of  penetration, 
istance  from  nose  tip  to 
le  from  target)  are  shown  in 
Case  1 is  seen  to  be  approxi- 
i.e.,  the  decelerations  are 
enetrator  weights  (or  W/A) . 
ned  in  detail  from  Figure  5, 


which  is  a plot  of 


penetration  depth 
nose  length 


where  the 
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subscripts  refer  to  the  case  number.  The  scaled  deceleration 
ratio  is  seen  to  oscillate  around  a value  of  1 for  about  the 
first  lialf  of  nose  embedment;  the  ratio  then  steadily  increases 
to  a value  of  1.06  at  a penetration  depth  of  1-1/3  nose  lengths 
(8  in.).  The  gradually  larger  force  seen  in  Case  2 is  due  to 
its  penetration  velocity  staying  higher,  relative  to  Case  1. 

At  8 in.  penetration  depth,  the  penetration  velocities  were 
191  ft/sec  for  Case  1 and  197  ft/sec  for  Case  2. 

A series  of  field  plots  at  penetration  depths  of  2,  4, 

6,  and  8 in.  were  obtained  for  each  case.  These  plots  show 
the  overall  projectile  and  target  configuration,  the  extent 
of  plastic  deformation  in  terms  of  the  generalized  plastic 
strain  (e^)  , the  velocity  vectors,  and  the  principal  stress  field. 
These  plots  are  presented  in  Figures  6-13  for  Case  1 and 
Figures  14-21  for  Case  2. 

Comparative  plots  of  the  loadings  on  the  penetrator, 
i.e.,  distributions  of  normal  stress,  axial  force,  and 
radial  force,  for  penetration  depths  of  2,  4,  6,  and  8 in. 
are  shown  in  Figures  22-25.  Finally,  a comparison  of  the 
loadings  for  Case  2 at  8 in.  depth  and  for  Case  1 at  the 
same  penetration  velocity  as  Case  2,  i.e.,  197  ft/sec,  is 
shown  in  Figure  26.  The  penetration  deptli  at  wliich  the 
penetration  velocity  for  Case  1 was  at  this  value  (197  ft/sec) 
was  3.85  in . 

REFERENCE 

1.  M.  H.  Wagner,  K.  N.  Kreyenhagen,  and  W.  S,  Goerke,  Numerical 
Analysis  of  DNA  Earth  Penetrator  Experiment  at  DRES,  DNA 
Report  3537F,  June  1975. 
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TABLE  1.  PENETRATOR  PROBLEM  CONDITIONS 


Case 


Case 


Pro j ec  t i le 

Target 

Impact 
Veloc i ty 

Weight  D 

iameter 

W/A 

Nose 

200  ft/sec 

64125  gm 
(141.4  lb) 

6 in . 

5 ps i 

1.25  CRH 

Sand 

200  ft/sec 

192375  gm 
(424.1  lb) 

6 in . 

15  ps  i 

1.25  CRH 

Sand 

TABLi:  2.  PROPHRTIIiS  OF  SAND  TAROF^T 


Density, 

Air  fraction, 

cl 

Cohes ion 

Unconfined  compressive  strength 

Unconfined  tensile  strength 

Mises  limit,  i'  v/sTT 

’ max  2 


1 .4256  gm/cm ’ 
89  Ib/ft^ 

4 0?, 

0 

0 

0 

2906  psi 


Initial  value  of: 


Bulk  modulus,  K 
’ 0 

3292 

psi 

Shear  modulus, 

2469 

ps  i 

Constrained  modulus,  M 

o 

6585 

ps  i 

Young's  modulus, 

5926 

psi 

Poisson's  ratio,  v 

o 

0.2 

Dilatational  velocity, 

585 

ft/sec 

Figure  2.  Initial  Computational  Grid. 
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Figure  3.  Penetrator  Deceleration  vs  Depth  of  Penetration. 
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Figure  4.  Time  Histories  of  Penetrator  Deceleration,  Velocity,  Depth  of 
Penetration,  and  Attached  Length. 
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Figure  9.  Velocity  Field  and  I’rlnclpal  Stress  Field  at  k In.  Penetration  Depth,  Case 
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Figure  10.  Grid  Configuration  and  Generalized  Plastic  Strain  Field  at  6 in.  Penetration  Depth, 
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Figure  12.  Grid  Configuration  and  Generalized  Plastic  Strain  Field  at 
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Figure  14.  Grid  Configuration  and  Generalized  Plastic  Strain  Field  at  2 in.  Penetration 
Depth,  Case  2. 
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Figure  15.  Velocity  Field  and  Principal  Stress  Field  at  2 In.  Penetration  Depth,  Case 


Figure  16.  Grid  Configuration  and  Generalized  I’l.istic  Strain  Field  at  4 in.  I’enet  r.it  ion  Depth,  Case 


Figure  17.  Velocity  Field  and  Principal  Stress  Field  at  4 in.  Penetration  Depth,  Case 
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Figure  18.  Grid  Configuration  and  Generalized  Plastic  Strain  Field  at 


Figure  19.  Velocity  Field  and  Principal  Stress  Field  at  6 in.  Penetration  Depth,  Case 
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Figure  21.  Velocity  Field  and  Principal  Stress  Field  at  8 in.  Penetration  Depth,  Case 


Loading  Distributions  at  2 in.  Penetration  Depth. 


Figure  23.  Penetrator  Loading  Distributions  at  A in.  Penetration 


Penetrator  Loading  Distributions  at  6 in.  Penetration  Depth. 


Penetrator  Loading  Distributions  at  8 in.  Penetration  Depth. 


In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 
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